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Abstract. The weak nonlinear Kerr interaction between single photons and intense 
laser fields has been recently proposed as a basis for distributed optics-based solutions 
to few-qubit applications in quantum communication and computation. Here, we 
analyze the above Kerr interaction by employing a continuous-time multi-mode model 
for the input /output fields to/from the nonlinear medium. In contrast to previous 
single-mode treatments of this problem, our analysis takes into account the full 
temporal content of the free-field input beams as well as the non-instantaneous response 
of the medium. The main implication of this model, in which the cross-Kerr phase 
shift on one input is proportional to the photon flux of the other input, is the existence 
of phase noise terms at the output. We show that these phase noise terms will preclude 
satisfactory performance of the parity gate proposed by Munro, Nemoto, and Spiller 
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1. Introduction 

Years of experimental and theoretical research on quantum information science have 
revealed the numerous difficulties that must be overcome to build a large-scale quantum 
computer. Such a computer could factor large numbers and search unstructured 
databases far more efficiently than any classical computer. To pave the road to quantum 
computer implementations, we first need to master few-qubit technology. One appealing 
approach to the latter problem uses optics-based configurations for quantum processors. 
Light can carry quantum information in the form of single-photon polarization states, 
or vacuum plus single-photon superpositions. Spontaneous parametric down-conversion 
can produce entangled photon pairs and heralded single photons. Beam splitters and 
wave plates can be used to convert between single-photon polarization and vacuum 
plus single-photon qubits. They can also accomplish arbitrary single-qubit rotations. 
To complete a universal gate set for quantum computation, all that is then needed 
is an appropriate two-qubit quantum gate, such as a controlled-NOT (cnot) gate. A 
universal gate set will also serve quantum communication systems, in that it can realize 
their required full Bell-state measurement (BSM) device. 

Realizing an all-optical cnot requires a nonlinear interaction between single 
photons, but traditional nonlinear materials offer only weak coupling between single 
photons. Linear optics quantum computing (2j [3] circumvents this difficulty by 
using photodetection and post-selection to provide a strong nonlinearity at the 
cost of a non-deterministic architecture that needs substantial ancillary resources. 
Electromagnetically-induced transparency [3] may afford a strong cross-Kerr effect at 
the single-photon level, and lead to a deterministic all-optical cnot implementation. 
Munro, Nemoto, and Spiller pQ, however, have made a rather different suggestion for an 
all-optical cnot. By successive weak cross-Kerr interactions between a strong coherent- 
state probe beam — that acts as a quantum communication bus — and a pair of single- 
photon qubit beams they realize a deterministic parity gate from which a cnot can be 
constructed. Their parity gate can also be used to realize a full BSM apparatus. Munro, 
Nemoto, and Spiller use a single-mode treatment of cross-phase modulation (XPM) to 
show that their parity gate achieves near-ideal performance. It has been known for some 
time, however, that a proper quantum theory for self-phase modulation (SPM) requires 
a continuous-time multi-mode theory [5j. Moreover, we have recently extended the 
continuous-time quantum theory of SPM to the case of XPM |6j. Using that theory we 
proved that the phase noises needed to ensure preservation of the free-field commutator 
brackets at the output of the XPM interaction precluded high-fidelity operation of the 
quantum phase gate proposed by Chuang and Yamamoto [7j, which relies on a strong 
nonlinear interaction between two single photons in a cross-Kerr medium. The question 
that we will address in the present paper is whether a similar fidelity degradation will 
occur in the parity gate from pQ when its performance is assessed in the continuous-time 
framework. 

There are two issues that necessitate employing a continuous-time model for the 
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cross-Kerr effect. First, the photonic qubits at the input to and output from a quantum 
gate will be flying qubits, i.e., free- field optical pulses. Single-mode treatments of the 
quantum gate that ignore the temporal behavior of the optical pulses have assumed — 
either explicitly or implicitly — that gate operation occurs within a high-finesse cavity. 
However, even then free-field inputs and outputs must be properly accounted for so that 
the fidelity of a multi-gate quantum circuit can be determined. In other words, if we put 
our optical gate in a black box, then its input-output relation is a transformation from 
a continuous-time free-field operator at its input to another continuous-time free-field 
operator at its output. 

The second issue that requires a multi-mode treatment of XPM, is the non- 
instantaneous nature of the cross-Kerr interaction. Optical fibers have a 5-10 fs response 
time for the SPM interaction [5], and a similar response time is expected for XPM. The 
response time of the nonlinear material plays a key role in its quantum behavior when 
single-photon pulses are involved. For example, we have shown [U] that if the XPM 
response time is much shorter than the time duration of a pair of single-photon input 
pulses, then essentially no XPM phase shift will occur on these pulses even when the 
peak nonlinear phase shift from a single photon is the ir rad that Chuang and Yamamoto 
require for their single-photon controlled-phase gate [TJ. This counter-intuitive behavior 
occurs because the 7r-rad phase shift only affects a very short — and randomly-distributed 
in time — portion of the single-photon pulse. Conversely, in the slow-response regime — 
wherein the optical pulses are much shorter than the XPM response time — we can 
adjust the timing of the input single-photon pulses so that one photon can induce an 
appreciable XPM phase shift on the other. 

For an XPM medium whose response function is causal and non-instantaneous, 
preservation of the free-field commutator relations for the output field operators imposes 
phase noise terms on these outputs. These phase noises will be seen, in this paper, to 
severely degrade the fidelity of the parity gate proposed by Munro, Nemoto, and Spiller. 
It will also be shown that the mean-squared phase noise is proportional to the response 
function's amplitude, implying that stronger nonlinearity is accompanied by increased 
phase noise. 



In this paper, we first summarize the model introduced in [6] for the XPM 
interaction, which applies in the absence of loss, dispersion, and SPM. Then, in SectionEJ 
we apply the slow-response version of this model to a simple gate in which a single- 
photon pulse induces a weak phase shift on a coherent-state probe beam. This gate 
is the building block for the Munro, Nemoto, and Spiller parity gate, which will be 
described and analyzed in Section HI That treatment constitutes the continuous-time 
fidelity analysis of their parity gate. Section [5] includes some numerical results, which 
quantify the parity gate's phase-noise induced fidelity loss. Section E] concludes the 
paper by discussing the applicability of our results to different scenarios. 
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2. Continuous-time cross-phase modulation 

In this section, we describe the continuous-time model for XPM, introduced in [6]. This 
model applies primarily to the Kerr interactions that occur in a length of a macroscopic 
material, e.g., an optical fiber. Hence, its input and output field operators are free 
fields. Whereas discrete modes are proper choices for modeling light that is confined in 
a cavity, a continuous-time formalism is more appropriate for free-field operators. In this 
formalism, a positive-frequency photon-units field operator associated with a +2;-going 
electric field, in a well-defined polarization, can be written as follows [HI E] 

Mt) = /^«Me-" (1) 



where a(u) is the annihilation operator associated with frequency u that satisfies 
[a(u), d>(u')] = 5(uj — to'). The limits of integration in ([T]) are from to oo. However, 
for the optical sources we shall consider, whose bandwidths are narrow in comparison 
to their common center frequency ujq ^> 0, we can extend the integral's lower limit to 
— oo. It then follows that 

[E(t),E\t')]=S(t-t') (2) 

is the free-field commutator. 

For our later calculations it will be useful to write our field operators in terms of a 
discrete set of basis functions: 

£(f) = X)Mi(*), (3) 

i 

where {<pi\ is a complete orthonormal set of functions satisfying 

fdt(f> i (t)<f>* j (t) = 6 ij and J E<t>* i (t)<f> i (lf) = 5(t-lf), (4) 

i 

and di = Jdt 4>*(t)E(t) is a discrete-mode annihilation operator, which satisfies [a*, aj] = 
5ij. Equation provides us with a prescription for converting a continuous-time field 
operator to a sum of discrete-mode operators, where each mode has a pulse shape 
orthogonal to that of the other modes. We will use this formalism frequently in 
forthcoming sections. Strictly speaking, such a discrete representation should be applied 
to a time-limited or band-limited field, but inasmuch as we will only need one mode per 
field for our study of slow- response XPM, no loss of generality is entailed by use of ([3]). 

Now, let us begin our consideration of the Kerr medium. We will assume that 
there is no loss, no dispersion, and no SPM in this medium. Furthermore, we denote 
the photon-units input field operators by Es(t) for the signal beam and Ep(t) for the 
probe beam, whose respective photon-units output field operators are E' s (t) and E' P (t); 
see figure [TJ Then, the input-output relationship for this material is given by [0] 

E' K {t) = e iiK ®e i(lK ®E K (t), for K = S,P, (5) 

where we have suppressed the group delay, and 

fi K (t) = k fdTh(t-T)El(r)Ej(r), for J, K G {S, P}, J ^ K (6) 
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E' P (t) = e l ^% lMt) E P (t) 
E' s (t) = e^V«W£ 5 (t) 



Figure 1. Modeling the cross-Kerr nonlinearity using input/output field operators. 
Here, £5 and £p are phase-noise operators, while (is and P-p are XPM operators that 
depend on the probe-beam and signal-beam photon-flux operators, respectively, as well 
as the medium's response function. 



gives the XPM phase shifts on the signal and the probe beams. Here, k is the XPM 
coupling coefficient, and h(t) is a causal response function, which has been normalized 
to satisfy Jdth(t) = 1. Equation © is in accord with our semiclassical understanding 
of the Kerr effect, in which we assume that the nonlinear operators responsible for the 
Kerr effect are proportional to the photon-flux operators. This is what we expect to hold 
for an ideal cross-Kerr medium. However, for the output-field operators to commute 
with each other, i.e., 

[Em,E' P (t')} = [E' s (t),E'Ut > )} = 0, (7) 

as is the case for their corresponding input-field operators, we need to include a pair of 
Langevin noise operators. These noise operators represent coupling to localized noise 
oscillators that typically represent the molecular vibrations in the medium, and they 
result in Hermitian phase-noise operators £s(i) and £p(£) i n the field-operator input- 
output relations ([5]). The coupling coefficient between the light and the corresponding 
reservoir mode at frequency f2 turns out to be proportional to ^JllJjX}, where 

= fdth(t)sm(ttt) (8) 

must be non-negative for Q > 0, so that all the damping coefficients are positive. 
Furthermore, equation (J7|) requires that [6] 

[is(t), £ P (u)} = iK[h{u -t)- h{t - u)]. (9) 

In thermal equilibrium, and £p(t) can be taken to be in zero-mean joint Gaussian 
states with the following symmetrized correlation function 

(£ K (t)£ K (u) + £ K (u)£ K (t)) = 

k [— Hi(n)coth[hn/(2k B T)]cos[n(t-u)}, ioiK = S,P } (10) 

where kpT is the thermal fluctuation energy. 

The above model for XPM has several interesting implications. First, from 
equation (jTUI) . it is seen that the phase noise variance at time t, f° r K = S,P, is 

nonzero even if T — >• 0. This can result from photon-phonon interactions in the medium 
[3]. Second, it can be easily shown that an instantaneous response function in this 
model cannot reproduce the well-known classical results for XPM [6] . The fact that the 
response function has a nonzero effective time duration A becomes important when we 
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are dealing with single-photon pulse inputs, because the nonlinear phase shift imparted 
by a single-photon pulse depends on its duration being greater or smaller than A [5]. 
Let us elaborate more on this issue by considering the special case of a single-photon 
signal input in the slow-response regime. 

A single-photon signal pulse can be represented by the following state 

\l) s = Jdt<f>{t)\l t ) 3 , (11) 

where the wave function (j)(t) satisfies Jdt \4>{t)\ 2 = 1. Here, \lt)s is a multi-mode state 
that represents a single photon at time t; it satisfies Es(t')\l t )s = S(t — t')\0)s, where 
lO)^ is the multi-mode vacuum state for the signal beam. It follows that |0(£)| 2 is the 
probability density function for observing this photon at time t. Equivalently, we can 
use the discrete formalism from equation (jHJ) by employing a complete orthonormal basis 
in which <j>i(t) = (j)(t). Using such a basis, we have that is a state for which a\ is 
in the number state |l) ai) while the other {a,} are in their vacuum states. The average 
XPM phase-shift factor introduced by this single photon on the probe beam is therefore 

( e *AH*)) = Jdr \<P(r)\ 2 e iKh{t - T \ (12) 
where we have used 19 



cxp 



J dt g(t)E\t)E(t) = N jexp j dt (e 9 ® - l)E\t)E(t) 



(13) 



with M{f{E\ E)} denoting the normally-ordered form of the operator f(E\ E). 

In the slow-response regime, |0(r)| 2 behaves like 8{r — t ), relative to h(t — r), 
for an appropriate to near the pulse's center. In other words, the signal pulse's time 
duration, r , is much shorter than A, the time duration of the response function. From 
equation (fl2l . we then obtain 

/ e *Ap(*)\ _ e inh(t-t ) ^ /]^\ 

Now, if the probe pulse shaped is <p(t — t'), then the phase shift induced by our single 
photon on the probe's pulse at t — to + 1' is Kh(t'). We can maximize this phase shift 
by choosing t! = t h , where h{t h ) = max t [h(t)). 

Given that A « 1 — 10 fs for optical fiber, it is impractical to work in its slow- 
response regime. However, this is the only regime in which a useful quantum interaction 
may be seen. In the fast-response regime — in which, r ^> A — only a A-duration 
portion of the probe pulse undergoes a nonlinear phase shift, and the time location of 
this region is randomly distributed over the entire probe pulse with probability density 
function | </>(£) | 2 [6]. This property of the fast-response regime precludes its being useful 
for optics-based quantum computation. Therefore, we will limit our subsequent work 

I Interferometric measurements are crucial to optics-based quantum computation. Thus we will assume 
that the probe and the signal beams have identical pulse shapes except for a timing offset t'. This 
assumption is not essential, however, for our current discussion of the slow-response XPM. All we 
really need is for the probe pulse's effective duration to be comparable to that of the signal pulse, and 
hence much shorter than that of the response function. 
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to the slow-response regime. If promising results are obtained, it will then behoove us 
to find a Kerr medium with a more useful response timj§|. 

3. Kerr nonlinearity between a single photon and a coherent state 

Consider a Kerr medium with XPM coupling constant k and a response function h(t), in 
which there is no loss, no dispersion, and no SPM. Suppose we illuminate this medium 
with a signal pulse and a probe pulse, with the former being in a superposition of the 
vacuum state and the single-photon state with pulse shape <p(t + th), viz., 

\^n)s = a\0} s + (3 Jdt(f>(t + t h )\l t ) s , (15) 

with | ck | 2 + \j3\ 2 = 1, and the latter being in the coherent state \ap<f)(t)), for which 
Ep(t)\ap<f)(t)) = ap<p(t)\ap<p(t)) . In the mode-decomposition formalism, we can think of 
the signal pulse as being in the superposition state \ips) — a |0) as + /3|l)a S ; where \n) as is 
the Fock state associated with the modal annihilation operator as = Jdt <f>*(t + th)Es(t). 
Similarly, the probe beam can be taken to be in a coherent state \ipp) = \ctp) ap 
associated with the modal annihilation operator hp = Jdt <p*{t)Ep{t). All other input 
modes are in their vacuum states. The time shift between the signal and probe pulses 
results in the signal's inducing the maximum nonlinear phase shift on the probe pulse 
in the slow-response regime. The goal of this section is to determine the output density 
operator for this simple gate, which is a building block for the parity gate that will be 
studied in the next section. 

In order to find the density operator, we first find the normally-ordered 
characteristic functional for the output field operators E' s {t) and E' p (t), i.e., 

XN((s(t),(p(t)) = /e/ dKs( * ) ^(') e / d *^W^ t We-/ d '^ w ^ ( ' ) e-/ d '^W^(*)\. (16) 



The exact evaluation of the above functional is, in general, a tedious task. However, for 
operation within the slow-response regime, it is sufficient to employ the discrete-mode 
picture, characterized by the a$ and ap modes, for the input field operators, in which 

E s (t) = <p(t + th)as + vacuum-state modes (17) 
Ep{t) = 4>(t)dp + vacuum-state modes (18) 

because the slow-response condition ensures that the output field operators will then 
obey 

E' s (t) = 0(t + t h )e^ s{ - th) e ifls{ - th) d s + vacuum-state modes (19) 
E' p {t) = (p{t)e^ pi0) e iflp{0) ap + vacuum-state modes. (20) 

The two implicit assumptions in the above equations are t = and th > t , where 
to is the center of |0(i)|, and th is the delay between the signal and the probe beams 

§ In the above analysis, we have assumed that there is no SPM in our nonlinear material. It turns 
out that even if the SPM effect is present in the medium, it can be suppressed by operating in the 
slow-response regime. This SPM suppression occurs when h(th) ^ ^spm(O), where /ispm(0 is the SPM 
response function for the probe beam. 
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at which the maximum nonlinearity will be induced. Now, it is easy to see that if 
JdtQ(t)(j)(t + t h ) = and fdt(p(t)(j)(t) = 0, then the averaging in equation §Tjo§ is 
over the vacuum modes, whose normally-ordered characteristic functions are unity. 
Hence, the only functions for which the value of XN{(s{t), (p{t)) is nontrivial are 
(s(t) = Cs<K^ + th) an d (p(t) = Cp<t>{t). An equivalent characteristic function can 
then be obtained using the discrete-mode operators 

a' s = e^-^e^-^as (21) 
a' P = e^ V^ (0) a P , (22) 



where 



M-h) = ^jdr h(-t h - t)E p (t)E p {t) 
*Kh(-t h ) f &t E p {t)E p {t) 



= Kh(—th)a P ap = 0. (23) 

and similarly, 

£ip(0) ~ 8a s ds, where 6 = nhith)- (24) 

In both equations ( 1231) and ( |24|) we have dropped number-operator terms for modes that 
are in their vacuum states. 

Equation (l23l) reflects the fact that the probe beam does not phase shift the signal 
pulse, because it lags the signal pulse in time and the cross-Kerr effect is causal. On 
the other hand, the signal pulse induces a phase shift 9 on the probe beam, as shown 
in equation ( 1241) . Our new characteristic function will then read 

Xn(CsAp) = (e^se^re-^'re-^'s) 

= | Q ,|2/ e -2iIm{apCpexp[i| J >(0)]}\ 

_|_ |/3| 2 (1 — |£ 5 | 2 )( e - 2i WapCpexpK P (0)]cxp[j6l]}^ 

— a* j3( s (e~ 2ilra ^ ap ^p exp ^ p<y0 ^ e^ s ^~ th ^) , (25) 

where the averaging is taken over the phase-noise terms is{~th) an d £p(0). In order 
to perform the averaging in equation ( 1251) . we will rewrite the phase-noise terms in 
a normally-ordered form by introducing an annihilation operator b, whose quadrature 
components are the Hermitian operators b\ = Re{6} = is(-th)/ V29 and 62 = Im{6} = 
£p(0)/ \p2B. Because [£s(— U), £p(0)] = i9, we see that b satisfies [b,¥] = 1. Moreover, 
because £s and £p are in thermal states, that will also be the case for b. Thus, with 
N = (tfb), we get 

o 2 = (iU-h)) = (S(o)> 

= (2N + 1)6/2 

= k [— Hi(tt) coth[m/(2k B T)]. (26) 
J 2n 
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In our new formalism, we can replace (f(¥)g(b)) with (f(P*)g(P))p, where denotes 
statistical averaging over f3 = fli+ifa, which is a classical, zero-mean, isotropic, complex- 
valued Gaussian random variable with variance N. Using this fact, we then obtain for 
arbitrary complex parameters A and r\ 

(e Aif " s( -* h) e^ p(0) ) = e- A ^ /2 (e Ai « s e^ p ) ?SiCp , (27) 

where £s and £p are independent, identically distributed, zero-mean, real-valued 
Gaussian variables with common variance a 2 . An interesting observation is that for 
nonzero values of A and 77, there exists a nontrivial phase shift e~ Xr]ie ^ 2 , even when the 
phase-noise variance a 2 approaches zero. This is a consequence of the non-commuting 
nature of the phase-noise operators £s{—th) an d £p(0). 

The above formula makes it possible to handle the phase-noise averaging in 
equation (|25|) . In particular, it is easy to verify that the following density operator 
corresponds to the characteristic function Xn((s,(p) 

p SP (a,(3) = ( M 2 |0) a , s (0| ® \a P e^)a' p (a P e^\ 

+ \f3\ 2 \l) a , s (l\ ® \a P e^ + V) a , p (a P e^ + V\ 

+ a*Pe^\l) a/s (0\ ® \a P e^ +e ^) a/p (a P e^ +e ^\ 

+ a/3*e-^|0) 4 (l| ® \a P e^ +e ^) a>p {a P e^ +e ^\ )^ p . (28) 

The final averaging over the classical variables will be applied later, when we calculate 
the fidelity of the parity gate. It is important to note that the density operator predicted 
by a single-mode XPM theory, as used in [I] does not coincide with the above density 
operator evaluated at zero phase-noise variance, i.e., at £s = £p = 0. Whereas the 
first two terms in equation ( 1281) also appear in the density operator associated with 
the single-mode treatment, this is not the case for the last two terms. The reason for 
this difference is the extra phase term that appeared in equation ( |2"T|) because of the 
non-commuting nature of the signal and probe phase-noise operators. 

4. Parity-gate fidelity analysis 

A parity gate accepts two input qubits, in the general form \ip in ) = Pq\11) ab + 
/3i\10)ab + /#2|01)ab + Ps\00)ab, and it provides a classical outcome that heralds 
whether the output is the even-parity state /3 |H)ab + Ps\W)ab or the odd-parity 
state Pi\10) ab + $2|01)ab- Munro, Nemoto, and Spiller [T] cascaded two single- 
photon/coherent-state XPM interactions to produce the all-optical parity gate shown in 
figure [21 from which an all-optical cnot can be constructed. In this scheme, Alice (A) 
and Bob (B) encode their qubits in the horizontal (H) and vertical (V) polarizations 
of their respective single-photon pulses, but only Alice's H polarization and Bob's V 
polarization interact with the coherent-state probeljj] A single photon on Alice's H mode 

|| Figure 2 implicitly assumes that the cross-Kerr effect is polarization independent, so that the same 
probe-beam polarization can undergo XPM with Alice's horizontal polarization and Bob's vertical 
polarization. This assumption entails no loss of generality, because a wave plate can be inserted in the 
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Figure 2. An optics-based distributed parity gate that uses a weak cross-Kerr 
nonlinearity. Alice and Bob encode their qubits (single photons) in the horizontal 
and vertical polarizations of single-photon pulses. The polarizing beam splitters 
(PBS) guide Alice's horizontal polarization and Bob's vertical polarization to a 
Kerr medium in which they interact with a coherent-state probe beam. Each such 
interaction can induce a phase shift 9 on the probe beam. The —9 phase shifter 
deterministically changes the probe's phase, and the final beam splitter (BS), with 
near-unity transmissivity 77, models a displacement operator that is needed prior to 
measurement with a photon-number resolving detector. 



will impart a weak phase shift 6 <C 1 to the probe beam, as will a single photon on Bob's 
V mode. Munro, Nemoto, and Spiller [I] used a single-mode treatment of XPM to show 
that it should be possible to distinguish between even and odd parity states with high 
fidelity. In this section we will use continuous-time XPM theory, in the slow-response 
regime, to provide a more accurate fidelity analysis for their parity gate. Our conclusion 
about the viability of the figure El gate will be rather different from theirs. 

Before delving into our continuous-time fidelity analysis, it is germane to reprise 
the single-mode description of the parity gate, as presented in pQ. Suppose that the 
bus (probe) beam is in a coherent state \ap)p, and Alice and Bob are initially in the 
state Po\HH)ab + @\\HV)ab + @2\VH) ab + Ps\VV)ab- Each of these terms induces a 
different phase shift on the coherent mode. For instance, when Alice and Bob are in the 
state \HV)abi they impose a 29 phase shift on the probe beam so that its state becomes 
\ape 2ld )p. After the post-XPM —6 phase shift and the -«p field displacement provided 
by injecting the coherent state —^Jf}/ (1 — rf) atp at the highly-transmitting (1 — 77 1) 
beam splitter of figure El the probe is left in the coherent state \cip(e td — l))p. Similar 
calculations for all other input terms leads to the following Alice-Bob-probe output 
state: 

l^out) = (Po\HH) A b + fc\VV}AB)\0}p 

+ Px\ HV} A B\ap(e td - 1)) P + fo\V H) abW^ - 1)) P 

w ((3o\HH}ab + P3\VV)ab)\0}p 

+ f3 l \HV) A B\i0ap)p + !3 2 \VH)AB\-i0a P )p, for B < 1. (29) 

Assume that the number-resolving detector has unity quantum efficiency, and that the 
coherent-state strength is such that exp(— 6 ,2 |o;p| 2 ) 1. Then, when no photons are 

probe-beam path between the two cross-Kerr interactions to permit the use of polarization-selective 
XPM in realizing the distributed parity gate. 
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detected, we will conclude that Alice and Bob are in the even-parity state 

(3 \HH) AB + f3 3 \VV) AB 

WOl AB = / = (oUJ 

V^oP + l&l 2 

because \{0\±i9ap)\ 2 = exp(— 9 2 \ap\ 2 ) <C 1. Likewise, when n > 1 photons are detected, 
we will conclude that Alice and Bob are in the joint state 

p{n\ipout) 



i>n)AB - I 

y'tr[p(n|Vw)(Vwp)p] 

Pi p(n\t9a P )p\HV) AB + ft P (n\ - i9a P ) P \VH) AB , , , 
, , tor t/Cl 

y / tr[p(n|Vw)(Vw|™)p] 
p 1 \HV) AB + (-l) n (3 2 \VH) AB 



. for n > 1 (31) 

where \n) P is the probe's n-photon Fock state. When n is an odd integer, we can apply 
a 7r-rad phase shift to make all the n > 1 observations result in the odd-parity output 



state (fc\HV)AB + P 2 \VH) AB )/ ] /\p 1 \ 2 + \(3 2 \ 2 . 

A key figure of merit for the distributed parity gate is the success-probability 
distribution, {P n }, i.e., the probability of being in the desired state \ip n ) AB when n 
probe photons have been detected. The total success probability is then given by 



^success = J2 P n=J2 |p(w|AB<^n|^out)| 2 
n=0 n=0 

= -Pcvcn + -Podd) (32) 

where [•] denotes Bloch-sphere averaging, 



P cvcn = P = |p(0U(^ |Vw>| 2 = |A)| 2 + l^l 2 = 1/2 (33) 
is the average success probability for the even-parity case, and 

oo j _ e -6 2 \ap\ 2 

Podd = E P n = E |p(nUB(^n|^out)| 2 = ~ (34) 

n=l n=l 1 

is the average success probability for the odd-parity case. 

The preceding single- mode treatment shows that 9\otp\ ^> 1 yields near-unity 
success probability, whereas 9\ap\ 1 leads to P S ucccss ~ 1/2 indicating that it is 
then impossible to distinguish between the odd- and even-parity states in \ip ont ). The 
single-mode model, however, does not account for the causality-induced phase noise 
that is intrinsic to continuous-time XPM theory. To evaluate the fidelity of the parity 
gate, for the slow-response-regime continuous-time theory, we employ the formalism 
introduced in the previous section. For simplicity, we assume Alice and Bob are in a 
normalized tensor-product state \ip A )\ip B ), where 

\rl>A) = a\0)A + P\l)A and \^ B ) = a'|0) B + P'\1)b, (35) 

with 

|0) x = \0}k and \1) K = J dt (j>(t + t h )\l t ) K for K = A,B. (36) 
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Here, \1)k, for K = A, B, represents the qubit that has nonlinear interaction with the 
probe beam, and, for simplicity, we have assumed that there is no propagation delay 
between A and B. Hence, according to figure El the following correspondence holds 
between the number of excitations and the polarizations of Alice's and Bob's qubits: 



(37) 



\0) a =\V)a,\1)a = \H) a and |0) B = \H) B , \1)b 

The probe beam is in the time-shifted coherent state \cxp<p(t))p that undergoes the 
maximum nonlinear phase shifts from Alice's and Bob's qubits. Then, the output density 
operator for Alice, Bob, and the probe can be obtained by using the density operator 
in equation (l2"Hj) twice; once for the Alice-probe interaction and once for the Bob-probe 
interaction. Note that the density operator for the state associated with field operators 
E' A {t) and E' P {t) is given exactly by equation ( 1281) . by replacing the S subscripts with 
A's. At that point in the figure El setup, the probe beam is in a superposition of four 
different coherent states, given by equation ( 1281) . For each of these coherent states we 
can employ equation (128"]) again to obtain its corresponding output density operator. 
The last step is to replace ape 1 " 1 with ap(e ! ' 7_s ' — 1) to incorporate the effects of the 
phase shifter and the displacement operation. The final density operator is then given by 
the following equation, in which the states are associated with the output annihilation 
operators d' A = fdt (jf{t + t h )E' A (t), a' B = J dt (jf(t + t h )E' B (t), and d" P = fdt <p*{t)E P {t): 



PAP B (a,(3,a',P') = (|a| 2 |0) AA (0|pgg + |/3| 2 | 1)^(1 ^ 



+ aW CA |l)AA(0|pS ) + a/5*e-^|0) AA (l|pgJ) 



(38) 



where 



.(00) 

Ppb 



= |«'| 2 |0)bb(0| <g> 
+ \(3'\ 2 \1} BB (1\® 
+ a'*/3 , e* B \l) BB (0\ 
+ a'(3'*e-^\0) BB (l 



AZp-8) 



ap{e 

a P (e lip - 1 
apie 



pp 



atp(e 



- r 



pp 

»(£p-0/2) 



ap(e 
- 1 



pp 



ap{e 



i(( P -e/2) 



.(ii) 
Ppb 



apie 



i(£ P -0/2) 



PP 



apie 



i(£p-0/2) 



(39) 



= |a'| 2 |0) BB (0| <g> a P {e 
+ \(3f\l) BB (l\^\a P (e^ + ^ 
+ a'*(3'e^\l) BB (0\® 
+ a'P'*e~^\0) BB (l\( 



pp 
1 



a P {e^ p - 1 



apie 



pp 

i(ip+e/2) _ l 



ap{e 



pp 



apie 



i(tp+&/2) 



r 



.(10) 

Ppb 



.(01) 

Ppb 



= \a'\ 2 \0) BB (0\® 
+ \P'\ 2 \1} BB (1\® 
+ a'*P'e^\l) BB (0\ 
+ a'p'*e-^ B \0) 



a.p\e" v " r — 1 
ap ( e ^p+sm _ ! 



apie 



1 



pp 



pp 



pp 



pp 



apie 



i{tiP+e/2) 



(40) 



apie 



i(€p-«/2) 



apie 



ap(e* sr — 1 



i BB\ 



(1| ® a P (e lip - 1) 



pp 



apie 



i£p 
i£p _ i 



(41) 
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and the average in equation (13811 is taken over the phase noise terms £4, £p, and 
£p, which are statistically independent zero-mean Gaussian random variables with 
respective variances o~\ = a 2 B = a 2 and a 2 P = 2a 2 . 

Using the above density operator, we can obtain the average success probabilities 
for the even- and odd-parity cases as follows 

P ev en = fo(a P ,0)[C + De^ /2 

' ' - ' ' (42) 



1 + 2o-p\a P \ 



■.(C + De- a p /2 ), when a P < 1 



and 



Podd = C[1- f (a P , 6)} + De"*' 2 [f (V2a P , 0) - f (a P , 0) 



C 



+ De~ 



9 2 \ap\ 2 /(l+2a 2 p \a P \ 2 ) 



1 + 2a 2 P \a P \ 



4/2 



1 



1 



where 



fo(ap,0) 



1 + Aa 2 P \a P \ 



-(0|a P (e^ +e) - 



1 + 2a 2 P \a P \ 



<7p + 0<l, (43) 



I"p(£p+0)| : 



£p 



Op + 



< 1, 



9 2 |ap| 2 /(i+2 ( j| ) |ap| 2 ) 



1 + 2ap|ap| 



and 



C = 



\a(3'\ 4 +\a'/3\ l 
\a/3'\ 2 + \a'l3\< 



0.383 and D = 



2\a(3a'(3' 



i\2 



\a(3'\ 2 + \a'(3\' 



(44) 



0.117. (45) 



In this analysis we have assumed that Alice's and Bob's initial states are independent 
and uniformly distributed over their respective Bloch spheres. 

The parity gate's total success probability, P succe ss = -Peven + -Podd, measures how 
often it performs the required parity-separation task. Ideally, P S ucccss ~ lj something 
that only occurs when both of its constituents, P eve n and P dd ; approach their maximum 
values of 1/2. From equation (l4"2i) . we see that ideal P ev en behavior occurs when 
<7p|ap| 2 <C 1 and op 1. This contrasts with the single- mode theory for P eve n, from 
equation (1331) . in which P ev en is always 1/2. From equation (14"3"1) we see that the D term 
in P odd is negative but it vanishes when either crp|ap| 2 1 or <7p|ap| 2 ^> 1. We also 
have that the C term in P Q dd is positive, achieving its maximum value of 0.38 when 
# 2 |ap| 2 ^> 1. This contrasts with the single- mode theory, from equation ( 1341) . in which 
Podd approaches 1/2 for # 2 |o;p| 2 ^> 1. The reason for this difference is likely what we 
cited in equations (|271) and (1281) regarding the extra phase terms that are due to the 
non-commuting phase-noise operators. Overall, based on our continuous-time theory 
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Log w (a 2 P ) 

Figure 3. Success probabilities P CVC n, -Podd, and P SUC cgss versus phase noise variance 
for \a P \ = 100 and 9 = 20mRad. 

for the cross-Kerr effect, it seems that for the optimum performance of the gate we have 
to satisfy the following three conditions simultaneously: 

9 2 \a P \ 2 > 1, a 2 P \a P \ 2 < 1, op < 1. (46) 

In the next section, we numerically evaluate P eV en and -Podd — first with arbitrary choices 
of ap and Q\otp\ and then with a two-pole response function — to investigate whether 
these constraints can be satisfied. 

5. Numerical results 

In figures [3] and HI we have plotted the success probabilities P ev en, -Podd, and P SU cccss 
versus ap and 0|ap|, respectively. Here, we have assumed that it is possible to change 
the phase variance ap without affecting 9, even though both are related to the cross-Kerr 
medium's response function. Figure [3] shows that P ev0 n is very sensitive to the phase 
noise, dropping significantly if crp|ap| 2 > 1. On the other hand, P dd only drops slightly 
for moderate values of ap, and it regains its maximum value at high values of a p. Its 
maximum, however, is bounded below 1/2 to the value of C — 0.38. In figure HI it can 
be seen that P ev en is independent of 9, whereas P G dd requires 9\ap\ > n/2 to achieve 
reasonable performance. That observation justifies the use of 9\a>p\ = 2 in figures [3] and 

m 

In figures [3] and HI we assumed that we could independently vary 9 and ap in 
order to get a sufficiently large value for 9 2 /ap. Such is not likely to be possible in our 
cross-Kerr effect model because both ap and 9 = Kh(th) are functions of h(t). In order 
to increase the maximum of h(t) while keeping its area constant, we need to make it 
narrower. That however increases the bandwidth of H(Q), which may increase a p. On 
the other hand, apocn and 9 oc k. Therefore, in order to satisfy both # 2 |ap| 2 > 1 and 
<7p|ap| 2 1 conditions, one needs a high ratio of 9 2 /ap, or equivalently, a high value 
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Q\a P I 

Figure 4. Success probabilities P eV on, Podd^ and P S ucccss versus 8\ap\ for aj, = 10 
and |a P | = 100. 



0.5 




Lo# 10 | a P | 

Figure 5. Success probabilities P C vcn, -Podd, and P S uccoss versus the probe parameter 
\ap\ using the two-pole response function for the medium at T = OK and 70 — > 0. 



for k. This is not a desired behavior for a system that was designed to operate in the 
weak-nonlinearity regime. 

To clarify the above points, let's consider a concrete example. Suppose h(t) is the 
two-pole response function associated with the two-pole frequency response 

^) = / dtMt)e « = ?f -jL_, (47) 

where Qq is the vibrational resonance frequency of the medium and 70 is a damping 
rate. This function is a good fit to the actual Raman-gain spectrum of silica fibers in 
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the under damped single resonance regime < 7o/2 < Q , [TO], in which 



h(t) = ^ /n2 W ° \ fort>0. (48) 

V^o " 7o/4 

Previous work using this response function in the continuous-time XPM model has 
shown [S] that the best system performance occurs when T = K and 70 — > 0, in which 
case, we have 

POO 

o- 2 p = (2k/tv) dQHAQ) = Ktt . (49) 
Jo 

For 70 — * 0, we have h(th) = from which it follows that ap = 9, i.e., higher 
nonlinearity comes at the price of more phase noise. 

Figure plots the success probabilities P e ven, -Podd, and -Psuccess versus \ap\ for the 
underdamped two-pole response at T = K with 9\otp\ = 2. It can be seen that optimum 
system performance is achieved at a low value of |«p|, for which P succeS s is only about 
0.49. This agrees with what we observed before. For high values of \ap\, we have that 
cr|,|ap| 2 = 9 1 op 1 2 = 2 1 op I also has a high value, which, in turn, degrades the even-parity 
performance. For low values of \ap\, P Q dd — > because fo(ap,9) — > 1 in equation ( |43l) . 



6. Discussion and conclusion 



The results shown in figure [5] demonstrate that our causal, non-instantaneous model 
for the cross-Kerr effect precludes the distributed parity gate's achieving high success- 
probability operation. One should also bear in mind that the fidelity we have derived 
is for the ideal case of no loss, no dispersion, and no SPM. Each of these effects can by 
itself significantly degrade system performance. Moreover, the above performance has 
been achieved under the slow-response conditions, which makes its practical application 
to fiber-based XPM highly questionable. However, there are still two other issues that 
should be addressed before coming to a definite conclusion about the parity gate's 
feasibility. First, within the range of validity of our model, we may still be able to 
do better if we choose an optimum response function for our medium. The material's 
normalized response function h(t) must be causal and satisfy the following conditions 

roo 

/ dth(t) = l and F^O) > for Q > 0. (50) 
Jo 

The above constraints define a convex set of functions. In this set, we are interested in 
finding the function h(t) that maximizes the success probability of the gate. An easier 
problem is to find the function h(t) that maximizes the ratio 9 /a p. Both these problems 
are analytically, as well as numerically, difficult. Even if we knew such a function, it 
would be difficult to find a material with the desired response function. Nevertheless, it 
would be interesting to find the ultimate possible performance of the parity gate using 
our continuous-time model for the cross-Kerr effect. 

The second issue is the applicability of our model to atomic systems, possibly 
under electromagnetically- induced transparency conditions pE], illuminated by single 
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photons. There are several proposals that use such systems to provide an effective 
cross-Kerr nonlinearity for single photons [TTJ [TJl EH E]. Our model may or may 
not be applicable to such scenarios. In our model, we translate what we expect to 
occur classically in a pure cross-Kerr medium into quantum field-operator language, and 
then we make it self-consistent by introducing phase noise operators with appropriate 
commutators. This is not necessarily what happens when a photon interacts with a 
single atom or a small ensemble of atoms. In these atom-interaction scenarios we expect 
to get some pulse-shape broadening, but how to relate this effect to our medium's 
response function is yet to be investigated, although the work reported in [T51 [TBI IT7] 
may be of value in this regard. Moreover, if the physical reason behind the phase noise 
and the non-instantaneous response function is molecular vibrations, it is logical to ask 
how much vibrational noise a single photon may encounter when it interacts with a 
small number of atoms. In other words, if we don't expect that atomic vibrations are at 
all significant at the single-photon level, do we need to worry about phase noise or not? 
These concerns prevent us from immediately generalizing our assessment of the parity 
gate to all its possible implementations. Nevertheless, a full field-operator treatment 
of the nonlinear behavior of atomic systems in response to single-photon pulses needs 
to be developed to properly answer the preceding questions, as the analysis we have 
presented for continuous-time XPM suggests that initial assessments of XPM-based 
two-qubit optical gates may be overly (perhaps even wildly) optimistic. 
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